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Quarkonium States at 5-Factories 

E. Robutti 
I. N.F.N. - Sezione di Genova, Italy 

from the BaBar Collaboration 

An overview is given on recent progress in the study of quarkonium spectroscopy at the _B-factories. 
In particular, an updated status report is presented of the long list of "charmonium-like" resonances 
newly discovered, whose assignment as true charmonium states is in most cases at least controversial. 
Also, new measurements on the decay properties of bottomonium states above open-_B production 
thresholds are shown. 



I. INTRODUCTION 



Much of the progress attained in recent years in the 
study of the quarkonium spectra is owed to the mea- 
surements performed at i3-factories. The impressive 
amount of data recorded by the BaBar. and Belle ex- 
periments has made it possible to study rare decay 
chains and to look for as yet undiscovered resonances 
in the charmonium and bottomonium mass regions. 
Results presented here are based on different sub- 
samples of the full datasets recorded up to now by 
the two experiments, corresponding to integrated lu- 
minosities of about 430 fb^^ (BaBar - final) and about 
850 fb^^ (Belle). Significant contributions also come 
from the analysis of the various data samples recorded 
by the CLEO detector. 



II. THE CHARMONIUM REGION 

A. Charmonium production at B-factories 

Charmonium states are copiously produced at 
B-factories in a variety of processes. 
B mesons decay to charmonium in about 3% of cases. 
Charmonium can also be produced in the continuum 
e'^e~ — > cc process: for the J/ipj for instance, the 
inclusive production cross section is about 2.5 pb in 
the energy region of the B-factories. Of particular 
interest here is the "double-charmonium production" 
process, first observed by Belle, where a J/^ is 
produced together and exclusively with another 
charmonium state. 

States with C — +1 are formed in 77 fusion, where 
two virtual photons are emitted by the colliding 
e"'"e~pair {e~^e~ — > e"'"e~7*7* e~^e~{cc)). About 
1 million rjc are produced in 100 fb~^-equivalent of 
data via this mechanism. 

Finally, states with J^^ = 1^^ arc formed in 
initial state radiation (ISR) events, where a photon 
is emitted by the incoming electron or positron 
(e+e" 7(cc)). The cross section for ISR produc- 
tion of the J/tp at the energy of the T{AS) peak is 
about 2 pb. 



B. The A:(3872) state 

The discovery of the narrow AT (3872) state by Belle 
in 2003 [l| has marked the start of a "new" spec- 
troscopy in the charmonium mass region. Since then 
the resonant structure has been observed in differ- 
ent production and decay channels by CDF , DO, 
BaBar and again Belle. Besides the discovery produc- 
tion mode B X{3872)K, the X{3872) has been 
observed in prompt pp production 0]; searches in 
prompt e~^e~ production Q and in e~^e~OT 77 forma- 
tion [J] have given negative results so far. In addition 
to the discovery decay mode Ar(3872) J/iIj tt+tt" 
(consistent with AT (3872) J/tp p°), the decay chan- 
nels X(3872) ^ J/7/> 7 d and AT (3872) -> D^lf-K^ 
@ have been seen. 

Detailed angular analyses 0, together with con- 
straints on quantum numbers from observed produc- 
tion and decay modes, favour the J^*^ = 1++ assign- 
ment, although J^'^ = 2~+ is not ruled out. 
The "natural" interpretation of the Ar(3872) as a char- 
monium state meets with considerable difficulties 0: 
several alternative models have been proposed in re- 
cent years. Given the proximity of its mass to the sum 
of the D° and D*° masses, the X(3872) might be an 
S'-wave D^-D*° bound state, or "molecule" A dif- 
ferent approach interprets the state as a "tetraquark" 
structure, or a bound state of a diquark-antidiquark 
pair [lo| . In this case a total of 4 nearly-degenerate 
states (2 neutral -I- 2 charged) should exist in the same 
mass region. Other non conventional interpretations, 
such as a ccg hybrid, are disfavoured, mainly because 
of the mass, which is inconsistent with current predic- 
tions. 

In their most recent updates of the J/ijj tt+tt^ 
analyses, both BaBar. and Belle have shown re- 
sults for the neutral, as well as for the charged, B 
decay mode [Tlj. Using samples corresponding to 
413 fb^^ and 605 fb~^ respectively, they reconstructed 
B+ -> J/i) 7r+7r-is:+[34j] and B" ^ J/^/> Tr+Tr-iC^ 
decays (Fig. [1]). Belle reports the first statistically 
significant observation (5.9cr) of the production mode 
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FIG. 1: The J/ip -k^tx invariant mass distributions in the X(3872) region for selected B decay events: a) BaBar, 
B+ J/ij TV+TT-K+; b) BaBar, B" J/^j Tv+Tv-Ri; c) BeUe, B+ J/^P tv+tt-K+; d) BeUe, B° J/i: 7r+7r'if§. 



B° X(3872)ifg; weaker evidence is also seen m 
BaBar data. 

The comparison between the charged and neutral 
channels is of particular interest: in the tetraquark 
model, two different X states are produced in the two 
decays, their masses differing by ~ 7 10 MeV/c^; 
moreover, in D^-D*^ molecular models the ratio of 
the neutral to the charged B decay rate may show de- 
parture from unity Table |I] reports results for the 
ratio of branching fractions and for the mass difference 
AM = Mxii — Mx-i^' the ratio appears significantly 
higher than 0.1, while the measured values of AM are 
consistent with no mass splitting. 



TABLE I: Results from latest B — > J/ip tt^tt K analyses 
from BaBar and Belle. 





B(B"^X(3872)K") 
B(B+^X(3872)K+) 


AM(MeV/c2) 


BaBar. 


0.41 ± 0.24 ± 0.05 


2.6 ± 1.6 ± 0.4 


Belle 


0.82 ± 0.22 ± 0.05 


0.18 ± 0.89 ± 0.26 



Belle's observation of the X(3872) D^D de- 



cay mode has been more recently corroborated by a 
similar study by BaBar fl3i|. In a sample correspond- 
ing to 347 fb^^, 5+ (-8°) candidates have been re- 
constructed in the -> D"D*°K+{K^) de- 
cay mode, where the D'^ arc reconstructed in the 
D'^ K^TT^ , K^tt^tt'^ , K^TT^TT^TT^ chamicls, and 
combined with tt*^ and 7 to give D*^ candidates. The 

n — *0 

distribution of the D D invariant mass is shown in 
Fig. [21 along with the tt" invariant mass distri- 

bution from Belle's analysis. 

The measured value for the mass of the resonance is 
consistent with Belle's measurement; however, it is 
significantly higher than the average of the measure- 
ments from X(3872) J/tp tt+tt" decays (Fig. [3]). 
This might indicate the presence of two different 
states, as foreseen in the tetraquark model [14], decay- 
ing predominantly to different final states. A simple 
alternative explanation has been proposed in [l^ , ex- 
ploiting the fact that, because of the proximity of the 
]jO£)*o threshold, the peak position in the m{D'^D* ) 
spectrum is sensitive to the orbital angular momen- 
tum of the final state. 
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FIG. 2: The D^D (tt") invariant mass distributions in the X(3872) region for selected B decay events: a) Belle: 
^ D°D"Tr°K+{K^) ; h)BABm:B+ {B°) D^D'^K+iKl). 
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FIG. 3: Comparison of mass measurements for the 
X(3872) state at B-factories. 



C. New states at 3940 MeV/c^ 

In the mass region around 3940 MeV/c^, Belle has 
reported the observation of three resonant structures, 
provisionally named X(3940), r(3940) and Z(3930). 
The X(3940) has been seen as a peak in the re- 
coiling mass spectrum against a reconstructed J/ip 
produced in "continuum" e+e~interactions; its decay 
to DD* has been subsequently identified [l^. The 
y(3940) appears as a threshold enhancement in the 
J/ip to invariant mass in reconstructed B — > J/^ luK 
decays Finally, the Z(3930) has been seen in 



the DD invariant mass spectrum in selected 77 de- 
cays p. 

Despite the proximity of the measured masses (see Ta- 
bleini, the structures are usually associated to three 
different states, mainly because of the different pro- 
duction and decay mechanism in which they have been 
observed. The mass and width values, the observed 
decay channel and a dedicated angular analysis make 
the Z(3930) a good candidate for the 2^P2 state of 
charmonium, the Xc2(2P). For the other two states, 
proposed charmonium assignments include the r]c{SS) 
and Xci(2-P) for the X(3940), and the Xci(2P) and 
Xco(2P) for the ^(3940) [|. However, these assign- 
ments lying on much weaker bases, different, uncon- 
ventional interpretations (such as tetraquarks) have 
been proposed in this case too. 

BaBar. has recently confirmed the observation of the 
y(3940) in a similar study of the B+ J/ip ujK+ 
and B° J/ip ojKg decay modes [l^, conducted on 
a sample corresponding to 348 fb~^. Combinations 
of 7r"'"7r~7r'^ are selected if their invariant mass lies in 
the to mass region, and coupled to reconstructed J/ip 
and K+ (K^) to form B+ {B°) candidates. The he- 
licity angle between the 7t~^ and the tt^ directions in 
the TT+TT" rest frame is then used to apply a weight- 
ing procedure in order to separate the lo contribution 
from the tt~^ '!t~ tt'^ combinatorial component: this es- 
tablishes that practically all selected B mesons decay 
to a real uj. It is thus possible to ignore the shape of 
the TT+TT^Tr" invariant mass distribution (which, in the 
m{J /ip 7r+7r~7r°) region of interest, is distorted due to 
the proximity of the production threshold) , and to ex- 
tract the number of B signal events in different bins 
of the J/il) u! invariant mass (Fig. S]). A fit to the 
resulting distribution with the sum of a resonant and 
a smooth background shape yields the values of mass 
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and width reported in Table [TTl both lower and sig- 
nificantly more precise than those obtained by Belle. 
The measured ratio of the neutral to charged B decay 
rate to YK is Ry = 0.27 
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FIG. 4: BABMi: the J/il) uo invariant mass spectrum for 
selected B+ J/ip luK+ (above) and B° -> J/ip uK^ 
(below) events. Data in the neutral channel have been 
normalized to take into account the missing Ki^ and Ks 

contributions. The superimposed curves show the 
result of the fit. 



D. The Y vector states 

Another set of unexpected structures in the mass re- 
gion above the open charm production threshold has 
been revealed by the study of ISR events where a, J/4> 
or a ip{2S) is produced in the final state. 
The first of these states has been observed by 
BaBar. [23| as a broad enhancement around 
4.26 GeV/c^ in the spectrum of the J/tp tt+tt^ 
invariant mass (Fig. [5^) in ISR events e+e^ — > 
J/ijj 7r"'"7r~7isR. The selection required no other 
charged particle detected except those reconstructing 
the final state; the detection of the ISR photon, pref- 
erentially emitted along the beam direction, was not 
required. Moreover, the total transverse momentum 
and the inferred recoiling mass were required to be 
consistent with the kinematics of the ISR process. The 
provisional name given to the resonance is F(4260). 
In a similar analysis performed on a sample of 
13.3 fb~^ of data collected at energies between the 
T(15') and T(45') peaks, CLEO has been able to 
confirm the observation of the F(4260) [2l|. An 
independent confirmation has come from an energy 
scan performed by CLEO-c in the region 3.77 GeV < 



< 4.26 GeV, where the exclusive reaction e+e^ 
J/ip 7r"'"7r~ has been reconstructed. In the same anal- 
ysis, a significant enhancement at y/s — 4.26 GeV has 
also been observed in the J/tp tt^tt^ final state, thus 
providing evidence for the y(4260) J/4> ttOtt" de- 
cay mode 22] . 

By repeating on the V'(2S') 7r"'"7r~ final state the same 
kind of analysis leading to the observation of the 
y(4260), BaBar. has identified another broad struc- 
ture just above the production threshold, at a mass 
around 4.32 GeV/c^ (Fig.Eb) [H. 
More recently. Belle has repeated the same studies 
with much higher statistics. On a data sample corre- 
sponding to 548 fb^^, ISR events with J/ip tt+tt" final 
state have been reconstructed and the corresponding 
invariant mass spectrum analysed [IJ. Beside the 
clear peak of the y(4260) (Fig. [S};), another broad 
enhancement is visible, which Belle fits with a reso- 
nant shape and names y(4008). On a larger sam- 
ple (673 fb^^), the study has been repeated for the 
■0(25') 7r+7r~ final state [23]: the analysis of the result- 
ing spectrum confirms with high significance the pres- 
ence of the y(4350) peak (Fig.EJi), and reveals yet an- 
other broad structure, at a mass around 4.66 GeV/c^. 
The mass and width of all 'F' states, as obtained by 
Belle from fits to Brcit-Wigner resonant shapes, are 
reported in Table [TTl 

All these states share several properties: they are vec- 
tor states with the photon quantum numbers, J^^ = 
1^^; they have large total widths; they have only 
been observed in charmonium decay modes. Not only 
they do not correspond to any known charmonium 
state, but were they all real and distinct resonances, 
they would clearly overpopulate that mass region. 
The most favourable charmonium assignment for the 
y(4260) would be the V(4S') [2^], but the i/;(4415) 
has already been proposed as a candidate; moreover, 
a preliminary upper limit of 7.6 found for the ratio of 
the DD to J/tp TT+TT^ decay rate appears far too low 
for this assignment. Alternative explanations include 
tetraquark objects [23] and ccg hybrids [28l ]. 



E. Overview 

Table [ill summarizes the main properties of the 
states recently discovered at i?-factories in the char- 
monium region (it also includes the F(4160), seen 
by Belle in double-charmonium production and not 
discussed in this work). Except for the Z(3930), 
which is likely to be the Xc2(2P) state of charmo- 
nium, it is quite hard to fit them into the spectrum 
of conventional cc mesons. Some of them might in 
fact represent the first clear indication of a new non- 
conventional hadron spectroscopy. 
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TABLE II: Summary of the properties of recently discovered states in the charmonium region. 



State 


Production mode 


Decay mode 


Experiments 


Af(MeV/c ) 


r(MeV/c ) 


tPC 


X(3872) 


B 
PP 


-> XK 

— > X [anything] 




Belle, BaBar, 
CDF, DO 


3871.2 ±0.5 


< 2.3 


1++ 

2 ^ 


X(3875) 


B 


^ XK 


D D 71,00 


Belle, BaBar. 


3875.2 ± 0.7 


3.0 ti 4 ± 0.9 




X(3940) 




-> Jli>x 


OD* 


Belle 


3943 ± 6 ± 6 


< 52 


tP + 


y (3940) 


B 


^ YK 


T J 1 

J /ifJUJ 


T) 11 T) A T) An 

Belle, BaBar. 


o 1 A -l-3.8 1 1 r\ 

3914.6 I3 4 ± 1.9 


34 Ig ±5 


rf + 


Z(3930) 


77 


^ Z 


OD 


Belle 


3929 ± 5 ± 2 


29 ± 10 ± 2 


2++ 


y(4008) 




Y 




Belle 


4008 ± 40 til'" 


226 ± 44 ± 87 


1"- 


Y(4260) 




Y 




BaBar,, Cleo, Belle 


4247 ± 12 til 


108 ± 19 ± 10 


1— 


y(4350) 


e"'"e~ 


-> Y 


ip{2S)TT+n-, 


BaBar, Belle 


4361 ± 9 ± 9 


74 ± 15 ± 10 


1"- 


F(4660) 




-> Y 




Belle 


4664 ± 11 ± 5 


48 ± 15 ± 3 


1"" 


F(4160) 


€^ e~ 


J/ipY 


O'D* 


Belle 


4433 ± 4 ± 2 


+18 +30 
-13 -13 






FIG. 5: Evidence for the vector 'Y' states in ISR events at B-factories. Spectra for final state invariant masses are 
shown: a) BaBar, miJ/ip tt+tt"); b) BaBar, m{i(){2S) tt+tt"); c) Belle, m{J/ip tt+tt"); d) Belle, m(V'(2S) tt+tt"). 



III. THE BOTTOMONIUM REGION 
A. Bottomonium spectroscopy at B-factories 



their energy being tuned to the peak of the T{AS) res- 
onance, which decays in almost 100% of cases to a BB 
pair. 



The B-factories are not usually considered as ideal 
facilities for the study of the bottomonium spectrum. 
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However, the huge amount of data accumulated by 
BaBar and Belle in several years of running allows for 
the search for rare, non-BB decay modes. Moreover, 
the ISR mechanism provides a way to access directly 
the states below the open-i? threshold, making for 
some of the largest data sets ever recorded in this re- 
gion. 

Another important option is given by the possibility 
of running the e^e^ colliders at energies different 
from the T(4S') mass. In fact, both BaBar. and Belle 
have recorded, in short periods of time, data samples 
at various energies in the bottomonium region which 
outrun all pre-existing sample in that range. They are 
summarized in Table ITTll [35j 



TABLE III: BaBar, and Belle data samples away from 
the T(45') peak. 





Energy 


f Cdt (fb-^) 




r{2S) peak 


15 


BaBar. 


T(35') peak 


30 




T{4S) - 11.3 GeV/c^ 


11.3 




T(35') peak 


3 


Belle 


T{5S) scan 


2 




T(55') peak 


21.7 



many of the modes studied. 

Figure [5^ shows the distribution of m{fi^fj,~) vs. 
AM^TT for selected /z"'"/i-~7r+7r~ events. After a fur- 
ther selection, the signal yield for each of the decay 
modes is obtained from fits to the AM^^r distributions 
for events within appropriate windows in m{£^£^). 
Extraction of the branching fractions provides the 
most precise determination to date of ;B(T(4S') 
T{1S,2S) n+n-) and 6(T(3S') ^ T(25) tt+tt-) and 
measurements of B{T{3S,2S) T(IS') tt+tt^) com- 
parable to current world averages. 
Figure shows the distribution of min^ n'^) vs. 
AM^ = m{e+e~n+n~n") - m(£+£-) - m(7r+7r-7r°) 
for selected ^+£^7r+7r^7r° events. In this case sig- 
nal yields are obtained from fits to the AM^ distri- 
butions for events within an appropriate window in 
m(7r'*'7r~7r°) (Fig. Significant yields (llcr and 

6.2a, respectively) are obtained in both the ji^ IJ,~ 
and e+e" samples, thus providing the first obser- 
vation of the T(4S') T(IS') ry decay mode. The 
resulting branching fraction, B{T{AS) T(IS') 77) = 
1.96 ± 0.06 ± 0.09, is about 2.4 times larger than the 
corresponding dipion rate. 



C. Anomalous T(55') decay rates 



B. Hadronic decays of the T(4S) 

Study of the hadronic transitions between bottomo- 
nium states provides a way to test QCD multipole 
expansion (QCDME) models, which are usually em- 
ployed to describe this kind of process. 
Observation of non-BB hadronic decays of the 
T(4S') has been first pubhshed by BaBar l29|], fol- 
lowing evidence for the T(4S') T(15) tt+tt^ 
decay mode previously reported by Belle; all re- 
sults have been confirmed by a more recent update 
from Belle (s^l- In these analyses, events with a 
muon pair and a charged pion pair tt+tt^ 
were selected. Then the distribution of TO(/i+/i^) 
vs. AM^TT = m(/U+/x~7r"'"7r~) — m(/x+/i~) was stud- 
ied in order to identify decay chains of the type 
T{mS) T{nS) tt+tt-, T(n5) ^ ii+ . In ad- 
dition to T(35) ^ T(15,25) tt+tt" and T(25) ^ 
T(IS') TT+TT~ events from ISR production, used as 
control samples, significant signals for T(4S') 
T[1S) TT+TT- and T(4S') -> T(25) tt+tt" were found, 
with branching fractions of the order of 10"^. 
In a recent update [Slj , BaBmi has used a data sample 
corresponding to 348 fb~^and added the study of the 
£~^£~r] final state, where the rj is reconstructed in the 
rj tt+tt^tt" channel. Due to the increased statistics, 
the e+e" sample now provides significant signals in 



(Note: this topic is discussed in much more detail 
in Ref. f33| in these Proceedings.) 
Another study of hadronic T(mS') T{nS) tt+tt" 
transitions has been performed by Belle, using the 
data sample recorded at ^/s ~ 10.87 GeV, near the 
peak of the T{5S). Fig. [7^ shows the distribution of 
m{fi'^fi~) vs. AMtttt for selected fj,~^ fi~ it~ events; 
the horizontal bands indicate the T(15'), T(2S') and 
T(35') signal regions (bottom to top) for TO(/x+/i"). 
Fig. [7)3 shows the AM^t^ projections for events in the 
T(15) and T(2S') band, where clear peaks corre- 
sponding to the T(35) ^ T(15) TT+TT", T(2S') ^ 
T(IS') TT+TT" and T(55) ^ T(IS') tt+tt", T(55) ~> 
T(25) TT+TT" transitions can be clearly seen. Simi- 
lar peaks are found, with a somewhat smaller signif- 
icance, in the T(5S') T{3S) tt+tt" and T{5S) 
T(15) K+K- modes. 

However, when partial widths are extracted from the 
data, they result in values more than two orders 
of magnitude larger than the corresponding ones for 
T(2S'), T(3S') and T(4S') decays. This observation 
prompts the question whether the observed enhance- 
ment at m(e+e") ^ 10.87 GeV/c^ really corresponds 
to the bottomonium state T(5S), or is rather sugges- 
tive of a "hidden beauty counterpart" of the y(4260) 
state observed in the charmonium region, which shows 
a surprinsingly large decay rate to J/ip tt+tt". 
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IV. SUMMARY 

The _B-factories have made a significant contribu- 
tion to our knowledge of the charmonium states (ob- 
servation of the rjc and possibly of the X2c{'2P), new 
measurements of resonance parameters); but most of 
all, in recent years, with the discovery of many unex- 
pected and often puzzling states in the charmonium 
mass region, may have unveiled a whole new type of 



spectroscopy. 

The quest for new discoveries is now being extended 
to the bottomonium sector, where thanks to the large 
amount of data accumulated and to dedicated run set- 
tings, the possibility of finding new exciting results 
appears concrete and is in fact already manifesting 
itself. 
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Charge-conjugate modes are implied here and in the 
following. 

A few weeks after this Conference, from an analy- 
sis of the data sample recorded at the T{3S) peak, 
BaBar. has reported the first observation of the ground 
state of bottomonium, the 77;, (15). 
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FIG. 6: Study of T{nS) hadronic decay modes at BaBar.: a)scatter plot of m{fi'^ jj.^) vs. AMtt-k for selected ^~-k^-k~ 
events: the boxes indicate the different signal regions; b) distribution of m(7r^7r~7r") vs. AM^ for selected ^~ -k'^ -k~ -k^ 
(red circles) and e"'"e~7r^7r~7r° (black crosses) events: the boxes indicate the signal regions in the two variables; c) 
projection of the AM,, distribution for jx^ events in the m(7r+7r~7r'') signal region in (b). 



- 10.8 

-SJ 

> 10.6 
a 10.4 
^ 10.2 



10 E* 

9.8 
9.6 
9.4 
9.2 
9 



(a) \x.\mii candidates 




"o 350 

5 300 

(3 250 

° 200 

2 150 
100 

.£ 50 

S 



? 80 



0.4 0.6 0.8 1 1.2 1.4 

AM = M(mx:m)-M((i|j,) (GeV/c^) 

(a) 


I (a) jj^jiji cs 

\ \ 


ndidates in the ' 

Y(2S) 


Y(3S) 


region 

Y(4S) 

. 1 


"Y(5S)" 

1 


0.4 0.6 0.8 


1 1.2 1.4 

AM = l\/l(n^iit:t)-M((in) (GeV/c^) 




(b) (i(i:t:T Ci 
Y(3S) 


ndidates in tlie Y(2S)^n*tr region 
Y(4S) :'Y(5S)" 

- ..i 



1 

AM : 



1.2 1.4 

M([iniiii)-M(|in) (GeV/c^) 



(b) 



FIG. 7: Study of T(nS') tt+tt" final states in Belle's T(55') 
sample: a) scatter plot of m{fj,^jj.~) vs. AMnn for se- 
lected ^'^ fj,~ Tv'^ n~ events: horizontal bands indicate the 
T(15'), T{2S) and T{3S) signal regions; b) projected 
distributions in AM-^tt for events in the T(15') (above) 
and T(2S) (below) bands in (a). 



